Genome assemblies of two species of porcelain crab, Petrolisthes cinctipes and Petrolisthes manimaculis (Anomura: Porcellanidae).

Crabs are a large subtaxon of the Arthropoda, the most diverse and species-rich metazoan group. Several outstanding questions remain regarding crab diversification, including about the genomic capacitors of physiological and morphological adaptation, that cannot be answered with available genomic resources. Physiologically and ecologically diverse Anomuran porcelain crabs offer a valuable model for investigating these questions and hence genomic resources of these crabs would be particularly useful. Here, we present the first two genome assemblies of congeneric and sympatric Anomuran porcelain crabs, Petrolisthes cinctipes and Petrolisthes manimaculis from different microhabitats. Pacific Biosciences high-fidelity sequencing led to genome assemblies of 1.5 and 0.9 Gb, with N50s of 706.7 and 218.9 Kb, respectively. Their assembly length difference can largely be attributed to the different levels of interspersed repeats in their assemblies: The larger genome of P. cinctipes has more repeats (1.12 Gb) than the smaller genome of P. manimaculis (0.54 Gb). For obtaining high-quality annotations of 44,543 and 40,315 protein-coding genes in P. cinctipes and P. manimaculis, respectively, we used RNA-seq as part of a larger annotation pipeline. Contrarily to the large-scale differences in repeat content, divergence levels between the two species as estimated from orthologous protein-coding genes are moderate. These two high-quality genome assemblies allow future studies to examine the role of environmental regulation of gene expression in the two focal species to better understand physiological response to climate change, and provide the foundation for studies in fine-scale genome evolution and diversification of crabs.

The Complete Mitochondrial Genome of the Hermit Crab Diogenes edwardsii (Anomura: Diogenidae) and Phylogenetic Relationships within Infraorder Anomura.

A complete mitochondrial genome (mitogenome) can provide important information for gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few mitogenomes of hermit crabs (superfamily Paguridae) in the infraorder Anomura have been reported. This study reports the first complete mitogenome of the hermit crab Diogenes edwardsii assembled using high-throughput sequencing. The mitogenome of Diogenes edwardsii is 19,858 bp in length and comprises 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. There are 28 and six genes observed on the heavy and light strands, respectively. The genome composition was highly A + T biased (72.16%), and exhibited a negative AT-skew (-0.110) and positive GC-skew (0.233). Phylogenetic analyses based on the nucleotide dataset of 16 Anomura species indicated that D. edwardsii was closest related to Clibanarius infraspinatus in the same family, Diogenidae. Positive selection analysis showed that two residues located in cox1 and cox2 were identified as positively selected sites with high BEB value (>95%), indicating that these two genes are under positive selection pressure. This is the first complete mitogenome of the genus Diogenes, and this finding helps us to represent a new genomic resource for hermit crab species and provide data for further evolutionary status of Diogenidae in Anomura.

Comet assay in Aegla platensis (Decapoda: Anomura) using a non-lethal hemolymph field sampling for in situ monitoring of freshwater genotoxicity.

This study aimed to apply the comet assay on Aegla platensis crabs as a suitable non-destructive approach for in situ monitoring of freshwater genotoxicity. Animals were captured during four sampling periods in a stream under minor anthropogenic impacts in Southern Brazil. Crabs were captured with a hand net, then the hemolymph samples were collected, and the animals were released into the stream after a 20-min recovery time. Hemolymph samples were transported to the laboratory and used to perform the alkaline comet assay. Results showed an intermediate level in the DNA damage index (range 107.3-165.0 arbitrary unit). No significant differences were observed among the different sampling periods. Hemolymph was successfully used as a non-lethal source of biological samples, and the comet assay using A. platensis proved to be a feasible approach for genotoxicity studies.

Mitochondrial genomes of the land hermit crab Coenobita clypeatus (Anomura: Paguroidea) and the mole crab Emerita talpoida (Anomura: Hippoidea) with insights into phylogenetic relationships in the Anomura (Crustacea: Decapoda).

The infraorder Anomura is a species-rich clade of decapod crustaceans recognized by its remarkable disparity in terms of morphology, anatomy, ecology, physiology, and behavior. This study assembled and characterized the complete mitochondrial genomes of two anomuran species, the hermit crab Coenobita clypeatus and the mole crab Emerita talpoida. The AT-rich mitochondrial genomes of C. clypeatus and E. talpoida are 16,469 bp and 15,810 bp long, respectively, and are composed of 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. A 1,390 bp and 553 bp long intergenic space is assumed to be the D-loop in C. clypeatus and E. talpoida, respectively. Mitochondrial synteny in C. clypeatus is identical to that reported in other congeneric hermit crabs while synteny in E. talpoida is identical to that described for the confamilial mole crab Stemonopa insignis. No major differences occur between the studied species and their respective congeneric / cofamilial species in terms of nucleotide composition and codon usage profiles of PCGs. Selective pressure analysis in PCGs, rarely conducted in anomuran crabs, indicate that all these mitochondrial PCGs experience purifying selection and that this purifying selection is stronger in some (i.e., cox family genes and cob) compared to other PCGs (e.g., atp8). Most of the tRNA genes exhibited a typical 'cloverleaf' secondary structure with few exceptions in the two studied species. In C. clypeatus, tRNA-Ser1 lacks the thymine pseudouracil cytosine (TΨC) loop while tRNA-Phe and tRNA-Tyr each exhibit a deletion of the dihydroxyuridine (DHU) loop but not the arm. In turn, in E. talpoida, tRNA-Phe and tRNA-Arg exhibit a deletion of the DHU loop but not the arm while tRNA-Ser1 lacks the TΨC arm. A phylogenomic analysis based on translated PCGs confirms the monophyly of the infraorder Anomura and retrieves most/all relationships at the superfamily and family level previously reported for anomurans. The analysis supports the monophyletic status of the families Albuneidae, Lithodidae, Coenobitidae, and Porcellanidae. In turn, the superfamily Paguroidea, and the families Paguridae and Diogenidae are polyphyletic.

Taxonomic status and phylogenetic relationship of Anomura (Crustacea: Decapoda) based on mitochondrial sequences and gene order rearrangements.

The complete mitochondrial genome (mitogenome) contributed crucial information, which could improve the better understanding of molecular phylogenetic analysis, evolution, and gene rearrangements. However, only a few Coenobitidae and Albuneidae mitogenomes have been described, and controversies about the Anomuran phylogeny remain. Here, we determined three Coenobitidae species (Coenobita perlatus, Coenobita rugosus, Coenobita clypeatus) and one Albuneidae species (Emerita talpoida). As the representative species of genera or families, and even the first species of the genus, supplementing these four species is essential to study the phylogenetic relationships within Anomura. Most of the four novel mitogenomes have the typical 37 genes, except for E. talpoida, which lacked trnF. Through combining these new data with 31 Anomuran mitogenomes from Genbank, different gene rearrangement patterns and internal phylogenetic relationships of Anomura were investigated. In our phylogeny analysis, all the Anomuran species were clustered into one group, and the polyphyly of Paguroidea was well supported. Moreover, various peculiar mitochondrial gene orders (MGOs) were summarized among Anomura and preliminary determined their rearrangement mechanisms through CREx. Twenty different MGOs were suggested in our research, and we were focused seven MGO patterns (Pattern A-J) in the process of discussing the gene rearrangement mechanism. Currently making full use of the large taxon sampling, our results provide valuable information on the evolutionary status of Anomuran species and are available for systematic rearrangement and phylogenetic analyses.

Different gene rearrangements of the genus Dardanus (Anomura: Diogenidae) and insights into the phylogeny of Paguroidea.

Complete mitochondrial genomes (mitogenomes) can provide useful information for phylogenetic relationships, gene rearrangement, and molecular evolution. In this study, the complete mitogenomes of two hermit crabs, Dardanus arrosor and Dardanus aspersus, were sequenced for the first time and compared with other published mitogenomes of Paguroidea. Each of the two mitogenomes contains an entire set of 37 genes and a putative control region, but they display different gene arrangements. The different arrangements of the two mitogenomes might be the result of transposition, reversal, and tandem duplication/random loss events from the ancestral pancrustacean pattern. Genome sequence similarity analysis reveals the gene rearrangement in 15 Paguroidea mitogenomes. After synteny analysis between the 15 Paguroidea mitogenomes, an obvious rearranged region is found in D. aspersus mitogenome. Across the 13 protein-coding genes (PCGs) tested, COI has the least and ND6 has the largest genetic distances among the 15 hermit crabs, indicating varied evolution rates of PCGs. In addition, the dN/dS ratio analysis shows that all PCGs are evolving under purifying selection. The phylogenetic analyses based on both gene order and sequence data present the monophyly of three families (Paguridae, Coenobitidae, and Pylochelidae) and the paraphyly of the family Diogenidae. Meanwhile, the phylogenetic tree based on the nucleotide sequences of 13 PCGs shows that two Dardanus species formed a sister group with five Coenobitidae species. These findings help to better understand the gene rearrangement and phylogeny of Paguroidea, as well as provide new insights into the usefulness of mitochondrial gene order as a phylogenetic marker.

Molecular phylogenetic analysis of the Paguristes tortugae Schmitt, 1933 complex and selected other Paguroidea (Crustacea: Decapoda: Anomura).

Morphological characters, as presently applied to describe members of the Paguristes tortugae Schmitt, 1933 species complex, appear to be of limited value in inferring phylogenetic relationships within the genus, and may have similarly misinformed understanding of relationships between members of this complex and those presently assigned to the related genera Areopaguristes Rahayu McLaughlin, 2010 and Pseudopaguristes McLaughlin, 2002. Previously undocumented observations of similarities and differences in color patterns among populations additionally suggest genetic divergences within some species, or alternatively seem to support phylogenetic groupings of some species. In the present study, a Maximum Likelihood (ML) phylogenetic analysis was undertaken based on the H3, 12S mtDNA, and 16S mtDNA sequences of 148 individuals, primarily representatives of paguroid species from the western Atlantic. This molecular analysis supported a polyphyletic Diogenidae Ortmann, 1892, although incomplete taxonomic sampling among the genera of Diogenidae limits the utility of this finding for resolving family level relationships. Several hypotheses regarding the evolutionary relationships among hermit crab genera were refuted by the Kishino-Hasegawa (KH). Shimodaira-Hasegawa (SH) and Approximately Unbiased (AU) tree topology tests, among them the hypothesis that Areopaguristes is monophyletic. A lack of support for the monophyly of Areopaguristes calls into question the phylogenetic validity of gill number for the differentiation of Paguristes, Areopaguristes, and Pseudopaguristes. The study was inconclusive with regard to the relationships among these three genera, but previously unknown diversity within both Paguristes and Areopaguristes was demonstrated. Existence of an undescribed species confounded under the name Paguristes tortugae Schmitt, 1933 was supported by genetics, morphology, and coloration. A second undescribed species with remarkable similarity to Areopaguristes hummi Wass, 1955 was discovered based on genetics and coloration.

Aegla buenoi n. sp. (Decapoda: Anomura): first record of aeglid crab from Cinzas River basin, Brazil.

A new species of freshwater anomuran crab, Aegla buenoi n. sp., is described. The new taxon was collected from two streams within the Cinzas River basin, Paran state, Brazil. We used morphological and molecular data (COI mtDNA) to distinguish the new species from its congeners. Aegla buenoi n. sp. is differentiated by morphological diagnostic features of the cephalothorax, chelipeds, second abdominal epimeron, and uropods. Molecular results confirm the separation of A. buenoi n. sp. from closely related species (A. castro Schmitt, 1942, A. lata Bond-Buckup Buckup, 1994, and A. jacutinga Marl Teixeira, 2020). Hence, our study increases the known diversity of aeglids and reports the first species of Aegla from the Cinzas River basin.

ß-elimination of hyaluronate by red king crab hyaluronidase.

Crustacean hyaluronidases are poorly understood both in terms of their enzymatic properties and in terms of their structural features. In this work, we show that the hepatopancreas homogenate of the red king crab has a hyaluronidase activity that is an order of magnitude higher than its commercial counterpart. Zymography revealed that the molecular weight of a protein with hyalorunidase activity is 40-50 kDa. Analysis of the hepatopancreas transcriptome and results of cloning and sequencing of cDNA revealed a hyaluronidase sequence with an expected molecular weight of 42.5 kDa. Further analysis showed that hyaluronat enzymatic cleavage follows the [Formula: see text]-elimination mechanism, which is well known for bacterial hyaluronidases. The results of ion-exchange chromatography showed that the final product of hyaluronate degradation is unsaturated tetrasaccharide. Thus, we identified a new hyaluronidase of higher eukaryotes, which is not integrated into the modern classification of hyaluronidases.

Chromosome-level genome assembly of Paralithodes platypus provides insights into evolution and adaptation of king crabs.

The blue king crab, Paralithodes platypus, which belongs to the family Lithodidae, is a commercially and ecologically important species. However, a high-quality reference genome for the king crab has not yet been reported. Here, we assembled the first chromosome-level blue king crab genome, which contains 104 chromosomes and an N50 length of 51.15 Mb. Furthermore, we determined that the large genome size can be attributed to the insertion of long interspersed nuclear elements and long tandem repeats. Genome assembly assessment showed that 96.54% of the assembled transcripts could be aligned to the assembled genome. Phylogenetic analysis showed the blue king crab to have a close relationship with the Eubrachyura crabs, from which it diverged 272.5 million years ago. Population history analyses indicated that the effective population of the blue king crab declined sharply and then gradually increased from the Cretaceous and Neogene periods, respectively. Furthermore, gene families related to developmental pathways, steroid and thyroid hormone synthesis, and inflammatory regulation were expanded in the genome, suggesting that these genes contributed substantially to the environmental adaptation and unique body plan evolution of the blue king crab. The high-quality reference genome reported here provides a solid molecular basis for further study of the blue king crab's development and environmental adaptation.

Transcriptome profiling and in silico detection of the antimicrobial peptides of red king crab Paralithodes camtschaticus.

Endogenous antimicrobial peptides (AMPs) are evolutionarily ancient factors of innate immunity, which are produced by all multicellular organisms and play a key role in their protection against infection. Red king crab (Paralithodes camtschaticus), also called Kamchatka crab, is widely distributed and the best known species of all king crabs belonging to the family Lithodidae. Despite their economic importance, the genetic resources of king crabs are scarcely known and no full-genome sequences are available to date. Therefore, analysis of the red king crab transcriptome and identification and characterization of its AMPs could potentially contribute to the development of novel antimicrobial drug candidates when antibiotic resistance has become a global health threat. In this study, we sequenced the P. camtschaticus transcriptomes from carapace, tail flap and leg tissues using an Illumina NGS platform. Libraries were systematically analyzed for gene expression profiles along with AMP prediction. By an in silico approach using public databases we defined 49 cDNAs encoding for AMP candidates belonging to diverse families and functional classes, including buforins, crustins, paralithocins, and ALFs (anti-lipopolysaccharide factors). We analyzed expression patterns of 27 AMP genes. The highest expression was found for Paralithocin 1 and Crustin 3, with more than 8,000 reads. Other paralithocins, ALFs, crustins and ubiquicidins were among medium expressed genes. This transcriptome data set and AMPs provide a solid baseline for further functional analysis in P. camtschaticus. Results from the current study contribute also to the future application of red king crab as a bio-resource in addition to its being a known seafood delicacy.

Novel gene rearrangement in the mitochondrial genome of Coenobita brevimanus (Anomura: Coenobitidae) and phylogenetic implications for Anomura.

The complete mitochondrial genomes (mitogenomes) can indicate phylogenetic relationships among organisms, as well as useful information about the process of molecular evolution and gene rearrangement mechanisms. However, knowledge on the complete mitogenome of Coenobitidae (Decapoda: Anomura) is quite scarce. Here, we describe in detail the complete mitogenome of Coenobita brevimanus, which is 16,393 bp in length, and contains 13 protein-coding genes, two ribosomal RNA, 22 transfer RNA genes, as well as a putative control region. The genome composition shows a moderate A + T bias (65.0%), and exhibited a negative AT-skew (-0.148) and a positive GC-skew (0.183). Five gene clusters (or genes) involving eleven tRNAs and two PCGs were found to have rearranged with respect to the pancrustacean ground pattern gene order. Duplication-random loss and recombination models were determined as most likely to explain the observed large-scale gene rearrangements. Phylogenetic analysis placed all Coenobitidae species into one clade. The polyphyly of Paguroidea was well supported, whereas the non-monophyly of Galatheoidea was inconsistence with previous findings on Anomura. Taken together, our results help to better understand gene rearrangement process and the evolutionary status of C. brevimanus and lay a foundation for further phylogenetic studies of Anomura.

Transcriptomic analysis reveals insights into deep-sea adaptations of the dominant species, Shinkaia crosnieri (Crustacea: Decapoda: Anomura), inhabiting both hydrothermal vents and cold seeps.

BACKGROUND: Hydrothermal vents and cold seeps are typical deep-sea chemosynthetically-driven ecosystems that allow high abundance of specialized macro-benthos. To gather knowledge about the genetic basis of adaptation to these extreme environments, species shared between different habitats, especially for the dominant species, are of particular interest. The galatheid squat lobster, Shinkaia crosnieri Baba and Williams, 1998, is one of the few dominant species inhabiting both deep-sea hydrothermal vents and cold seeps. In this study, we performed transcriptome analyses of S. crosnieri collected from the Iheya North hydrothermal vent (HV) and a cold seep in the South China Sea (CS) to provide insights into how this species has evolved to thrive in different deep-sea chemosynthetic ecosystems. RESULTS: We analyzed 5347 orthologs between HV and CS to identify genes under positive selection through the maximum likelihood approach. A total of 82 genes were identified to be positively selected and covered diverse functional categories, potentially indicating their importance for S. crosnieri to cope with environmental heterogeneity between deep-sea vents and seeps. Among 39,806 annotated unigenes, a large number of differentially expressed genes (DEGs) were identified between HV and CS, including 339 and 206 genes significantly up-regulated in HV and CS, respectively. Most of the DEGs associated with stress response and immunity were up-regulated in HV, possibly allowing S. crosnieri to increase its capability to manage more environmental stresses in the hydrothermal vents. CONCLUSIONS: We provide the first comprehensive transcriptomic resource for the deep-sea squat lobster, S. crosnieri, inhabiting both hydrothermal vents and cold seeps. A number of stress response and immune-related genes were positively selected and/or differentially expressed, potentially indicating their important roles for S. crosnieri to thrive in both deep-sea vents and cold seeps. Our results indicated that genetic adaptation of S. crosnieri to different deep-sea chemosynthetic environments might be mediated by adaptive evolution of functional genes related to stress response and immunity, and alterations in their gene expression that lead to different stress resistance. However, further work is required to test these proposed hypotheses. All results can constitute important baseline data for further studies towards elucidating the adaptive mechanisms in deep-sea crustaceans.

Development and characterization of 20 polymorphic microsatellite loci in the deep sea squat lobster, Munida isos Ahyong and Poore, 2004 and cross-amplification in two congeneric species.

Munida isos is a deep sea squat lobster species that is widely distributed across the New Zealand and east Australian region, and is often associated with deep sea vulnerable marine ecosystems. To investigate its population genetic structure and patterns of regional connectivity, microsatellite loci were developed for M. isos from two genomic libraries using the Illumina HiSeq 2500 sequencing platform. Twenty-six loci amplified consistently in M. isos from the Tasman Sea, among which 20 were polymorphic and selectively neutral. Evidence of null alleles was observed at eight loci. Most loci exhibited moderate to high levels of polymorphism, with an average polymorphic information content value of 0.482. The mean number of alleles per locus was 7.45, with a mean expected heterozygosity of 0.520. Thirteen loci exhibited significant deviation from Hardy-Weinberg equilibrium, while only one locus pair was in linkage disequilibrium after false discovery rate correction for multiple testing (P < 0.05). Cross-species amplification tests revealed that the transferability of 14 loci (70%) was positive for the two congeners M. endeavourae and M. gracilis. The accessibility to new polymorphic microsatellite loci will facilitate population genetic studies and aid in developing conservation and management strategies for vulnerable marine ecosystems.

Large-scale mitochondrial gene rearrangements in the hermit crab Pagurus nigrofascia and phylogenetic analysis of the Anomura.

Complete mitochondrial genome (mitogenome) provides important information for better understanding of gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few Paguridae mitogenomes have been reported. Herein, we described the complete mitogenome of hermit crab Pagurus nigrofascia. The total length was 15,423 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNA genes, as well as an AT-rich region. The genome composition was highly A + T biased (71.4%), and exhibited a negative AT-skew (-0.006) and GC-skew (-0.138). Eight tRNA genes, two PCGs and an AT-rich region found to be rearranged with respect to the pancrustacean ground pattern gene order. Duplication-random loss and recombination model were adopted to explain the large-scale gene rearrangement events. Two phylogenetic trees of Anomura involving 12 families were constructed. The results showed that all Paguridae species were clustered into one clade except Pagurus longicarpus, which for the first time imposed raises doubt about the morphological taxonomy of this species. Furthermore, the present study found that higher- level phylogenetic relationships within Anomura were controversial, compared with the previous studies. Our results help to better understand gene rearrangements and the evolutionary status of P. nigrofascia and lay foundation for further phylogenetic study of Anomura.

MicroCT imaging applied to description of a new species of Pagurus Fabricius, 1775 (Crustacea: Decapoda: Anomura: Paguridae), with selection of three-dimensional type data.

A new species of hermit crab, Pagurus fraserorum n. sp. (family Paguridae) is described from rocky subtidal reefs off KwaZulu-Natal, South Africa, and illustrated using both conventional drawings and colour photographs, and via three-dimensional (3D) X-ray micro-computed tomography (µCT). Because of the limitation µCT has in detecting very fine and soft structures, a novel approach of manually drawing setation and spinulation onto the two-dimensional images of the 3D visualizations was used to illustrate the pereopods. In addition, an interactive figure and rotation movie clips in the supplement section complement the species description, and the 3D raw data of the 3D type data are downloadable from the Gigascience Database repository. The new species is the sixth species of Pagurus Fabricius, 1775 reported from South Africa and is closely allied to the Indo-Pacific P. boriaustraliensis Morgen, 1990 and P. pitagsaleei McLaughlin, 2002, from which it differs by its shorter ocular peduncles, by the armature of the carpus of the right cheliped, and also in colouration. This study presents the first description of a hermit crab in which a majority of taxonomic details are illustrated through 3D volume-rendered illustrations. In addition, colour photographs and COI molecular barcodes are provided, and the latter compared to COI sequences of specimens from Western Australia previously identified as P. boriaustraliensis and of specimens of P. pitagsaleei from Taiwan, as well as to three additional South African members of the genus. The South African taxon was confirmed to be genetically distinct from all species tested.

Cloning and characterization of serpin from red king crab Paralithodes camtschaticus.

Serpins are a family of serine protease inhibitors that are involved in numerous physiological processes and are known to regulate innate immunity pathways. To advance our understanding of their role in P. camtschaticus, a commercially significant species, we cloned and characterized a serpin from this species, designated serpin PC, that has anticoagulant and anticomplement effects on human blood. We found that serpin PC is a secreted protein with a typical serpin-like primary structure that is similar to other known crustacean serpins. Recombinant serpin PC was found to have inhibitory activity against R/K-specific bovine cationic trypsin. The reaction proceeds through the formation of a stable covalent complex of peptidase with P1 residue R383 of serpin PC. This interaction is characterized by a relatively high overall inhibition constant kass=(2.3 ±â€¯0.7) × 106 M-1s-1 and an SI of 4.7 ±â€¯0.8. Protein localization by western blotting showed that serpin PC is present in the muscles and, to a lesser extent, the heart, whereas it is transcribed predominantly in hemocytes and the heart. Through peptidase activity profiling of hemocytes and plasma, we found that serpin PC inhibits at least two R/K-specific activities and showed that it inhibits phenoloxidase (PO) activity induction in hemocytes.

ORDER within the chaos: Insights into phylogenetic relationships within the Anomura (Crustacea: Decapoda) from mitochondrial sequences and gene order rearrangements.

The infraorder Anomura consists of a morphologically and ecologically heterogeneous group of decapod crustaceans, and has attracted interest from taxonomists for decades attempting to find some order out of the seemingly chaotic diversity within the group. Species-level diversity within the Anomura runs the gamut from the "hairy" spindly-legged yeti crab found in deep-sea hydrothermal vent environments to the largest known terrestrial invertebrate, the robust coconut or robber crab. Owing to a well-developed capacity for parallel evolution, as evidenced by the occurrence of multiple independent carcinization events, Anomura has long tested the patience and skill of both taxonomists attempting to find order, and phylogeneticists trying to establish stable hypotheses of evolutionary inter-relationships. In this study, we performed genome skimming to recover the mitogenome sequences of 12 anomuran species including the world's largest extant invertebrate, the robber crab (Birgus latro), thereby over doubling these resources for this group, together with 8 new brachyuran mitogenomes. Maximum-likelihood (ML) and Bayesian-inferred (BI) phylogenetic reconstructions based on amino acid sequences from mitogenome protein-coding genes provided strong support for the monophyly of the Anomura and Brachyura and their sister relationship, consistent with previous studies. The majority of relationships within families were supported and were largely consistent with current taxonomic classifications, whereas many relationships at higher taxonomic levels were unresolved. Nevertheless, we have strong support for a polyphyletic Paguroidea and recovered a well-supported clade of a subset of paguroids (Diogenidae + Coenobitidae) basal to all other anomurans, though this requires further testing with greater taxonomic sampling. We also introduce a new feature to the MitoPhAST bioinformatics pipeline (https://github.com/mht85/MitoPhAST) that enables the extraction of mitochondrial gene order (MGO) information directly from GenBank files and clusters groups based on common MGOs. Using this tool, we compared MGOs across the Anomura and Brachyura, identifying Anomura as a taxonomic "hot spot" with high variability in MGOs among congeneric species from multiple families while noting the broad association of highly-rearranged MGOs with several anomuran lineages inhabiting extreme niches. We also demonstrate the value of MGOs as a source of novel synapomorphies for independently reinforcing tree-based relationships and for shedding light on relationships among challenging groups such as the Aegloidea and Lomisoidea that were unresolved in phylogenetic reconstructions. Overall, this study contributes a substantial amount of new genetic material for Anomura and attempts to further resolve anomuran evolutionary relationships where possible based on a combination of sequence and MGO information. The new feature in MitoPhAST adds to the relatively limited number of bioinformatics tools available for MGO analyses, which can be utilized widely across animal groups.

Discrete phenotypes are not underpinned by genome-wide genetic differentiation in the squat lobster Munida gregaria (Crustacea: Decapoda: Munididae): a multi-marker study covering the Patagonian shelf.

BACKGROUND: DNA barcoding has demonstrated that many discrete phenotypes are in fact genetically distinct (pseudo)cryptic species. Genetically identical, isogenic individuals, however, can also express similarly different phenotypes in response to a trigger condition, e.g. in the environment. This alternative explanation to cryptic speciation often remains untested because it requires considerable effort to reject the hypothesis that the observed underlying genetic homogeneity of the different phenotypes may be trivially caused by too slowly evolving molecular markers. The widespread squat lobster Munida gregaria comprises two discrete ecotypes, gregaria s. str. and subrugosa, which were long regarded as different species due to marked differences in morphological, ecological and behavioral traits. We studied the morphometry and genetics of M. gregaria s. l. and tested (1) whether the phenotypic differences remain stable after continental-scale sampling and inclusion of different life stages, (2) and whether each phenotype is underpinned by a specific genotype. RESULTS: A total number of 219 gregaria s. str. and subrugosa individuals from 25 stations encompassing almost entire range in South America were included in morphological and genetic analyses using nine unlinked hypervariable microsatellites and new COI sequences. Results from the PCA and using discriminant functions demonstrated that the morphology of the two forms remains discrete. The mitochondrial data showed a shallow, star-like haplotype network and complete overlap of genetic distances within and among ecotypes. Coalescent-based species delimitation methods, PTP and GMYC, coherently suggested that haplotypes of both ecotypes forms a single species. Although all microsatellite markers possess sufficient genetic variation, AMOVA, PCoA and Bayesian clustering approaches revealed no genetic clusters corresponding to ecotypes or geographic units across the entire South-American distribution. No evidence of isolation-by-distance could be detected for this species in South America. CONCLUSIONS: Despite their pronounced bimodal morphologies and different lifestyles, the gregaria s. str. and subrugosa ecotypes form a single, dimorphic species M. gregaria s. l.. Based on adequate geographic coverage and multiple independent polymorphic loci, there is no indication that each phenotype may have a unique genetic basis, leaving phenotypic plasticity or localized genomic islands of speciation as possible explanations.

Morphological and molecular data reveal the cryptic diversity among populations of Aegla paulensis (Decapoda, Anomura, Aeglidae), with descriptions of four new species and comments on dispersal routes and conservation status.

The taxonomy of the Brazilian aeglid species Aegla paulensis Schmitt, 1942 from two disjunct hydrographic basins is revised using morphological and molecular data. Results show that six disjunctive populations of Aegla paulensis form a species complex. Aegla paulensis sensu stricto is redescribed and Aegla rosanae Campos Jr., 1998 is revalidated. The four remaining populations previously assigned to Aegla paulensis are now recognized as different species, namely Aegla  vanini n. sp., Aegla japi n. sp., Aegla jaragua n. sp. and Aegla jundiai n. sp. All new species are described and illustrated and are well supported by both morphological and molecular data. Aegla lancinhas Bond-Buckup & Buckup in Santos et al., 2015, which until recently was confounded with Aegla paulensis s. str., is supported as a valid species. A key to all members of the A. paulensis species complex is provided, and their phylogenetic and biogeographic relationships to other closely related species is discussed.